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Abstract 

We calculate the 0{a) improvement coefficient csw in the Sheikholeslami-Wohlert quark action 
for various improved gauge actions with six-link loops. We employ a conventional perturbation 
theory introducing the fictitious gluon mass to regularize the infrared divergence. Our results for 
some improved gauge actions are in agreement with those previously obtained with the Schrodinger 
functional method. 
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I. INTRODUCTION 



Recently the CP-PACS collaboration shows by a large scale of simulation that the hadron 
spectra in the quenched approximation systematically deviate from the experimentally ob- 
served ones both in the meson and the baryon sectors l| . It is now obvious that the next step 
is to incorporate the effects of dynamical quarks to reproduce the correct hadron spectra. 
With the current computational resources, however, unquenched QCD simulations are often 
restricted on lattices with the lattice spacing coarser than 0.1 fm while keeping the physical 
volume larger than 2 fm. 

A practical way to reduce the scaling violation effects is to employ the improved quark 
and gauge actions. For the quark part the 0(a) improved action proposed by Sheikholeslami 
and Wohlert 3] is now widely used. This action requires only one new term called a clover 
term. Although from a theoretical point of view the plaquette gauge action is already 0(a) 
improved, a comparative numerical study employing the various quark and gauge actions 
shows that the renormalization group (RG) improved gauge action reduces non-negligible 
0(a 2 ) errors[3j]. Moreover, JLQCD collaboration has recently reported that the first order 
phase transition observed in the three flavor QCD simulation with the 0(a) improved quark 
action and the plaquette gauge action, which is considered to be a lattice artifact, disappears 
once the gauge action is replaced by the RG improved one|^|. Thus the improvement of the 
gauge action is mandatory for the three flavor QCD simulation at the currently accessible 
lattice spacing. 

In this paper we determine the clover coefficient csw m the massless SW quark action 
up to one-loop order for various improved gauge actions including the DBW2 action^. 
Preparing for new improved gauge actions yet to come, we parameterize the value of csw 
as a function of the improvement coefficient of gauge action for later convenience. Another 
important purpose of the present calculation is to check the validity of the conventional 
perturbative method for the determination of the massless clover coefficient csw- Although 
previous calculations of csw are done by the twisted antiperiodic boundary conditions [6] 
or the Schrodinger functional method j''], we instead employ the conventional perturbation 
theory with the use of the fictitious gluon mass to regularize the infrared divergence, which 
has been applied successfully for the calculation of the renormalization constants and the 
improvement coefficients for the bilinear quark operators [8]. This method can be easily 
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implemented, within the standard knowledge of perturbation theory. Our results for some 
improved gauge actions are in agreement with those previously obtained with the Schrodinger 
functional method, which assures the validity of our conventional perturbative method. We 
are now extending this calculation of csw to the case of the heavy quark formulation proposed 
by the authors^], where the conventional perturbative method is much easier to handle the 
massive quarks than the Schrodinger functional method. 

This paper is organized as follows. In Sec. II we introduce the improved quark and 
gauge actions and their Feynman rules relevant for the present calculation. In Sec. Ill 
we determine the clover coefficient csw up to one-loop level from the on-shell quark-quark 
scattering amplitude. The result of csw is parametrized as a function of the improvement 
coefficient of the gauge action. Our conclusions are summarized in Sec. IV. 

The physical quantities are expressed in lattice units and the lattice spacing a is sup- 
pressed unless necessary. We take SU(iV c ) gauge group with the gauge coupling constant 
9- 

II. ACTION AND FEYNMAN RULES 

For the quark action we consider the 0(a)-improved quark action^]: 

n /x n 

f _ 

-cswY.Y. i 9j' l f J ^^F flu {n)tp n , (1) 

where we define the Euclidean gamma matrices in terms of the Minkowski matrices in the 
Bjorken-Drell convention: jj = -i>y J BD (J = 1, 2, 3), j 4 = j BD , 75 = j BD and = \[l^ 7,,]. 
The field strength F^ v in the clover term is given by 
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The weak coupling perturbation theory is developed by writing the link variable in terms of 
the gauge potential 

U nifl = exp [igaT A A A (n + \p)), (7) 

where T A (A = 1, . . . , N% — 1) is a generator of color SU(iV c ). 

The quark propagator is obtained by inverting Wilson Dirac operator in eq.Q, 
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. (k Xa + 1k 2il + k 3 A (q v -p u \ fhu- k u 
+o MT sm I — - ^- I cos ( ) cos 



where /abc the structure constant of SU(iV c ) gauge group. The first three vertices originate 
from the Wilson quark action and the last three from the clover term. The momentum 
assignments for the vertices are depicted in Fig. ^ 

For the gauge action we consider the following general form including the standard pla- 
quette term and six-link loop terms: 

S % = ~2 I C ° 51 tlU Pl + C l tlU rt g + C 2 ^2 tlU chr + C 3 ^ tvllplg > (15) 

9 I plaquette rectangle chair parallelogram J 



with the normalization condition 



c + 8ci + 16c 2 + 8c 3 = 1, 



(16) 



where six-link loops are composed of a 1 x 2 rectangle, a bent 1x2 rectangle (chair) 
and a three-dimensional parallelogram. In this paper we consider the following choices: 
O(Plaquette), a = -1/12, c 2 = c 3 = O(Symanzik) h], Q ci = -0.331, 



Cl 



C-2 



C 3 



13, ci = -0.252, c 2 + c 3 



c 2 = c 3 = O(Iwasaki), c\ = — 0.27, c 2 + c 3 = — 0.04(Iwasaki') 

-0.17(Wilson)Q and c x = -1.40686, c 2 = c 3 = 0(DBW2)[5]. The last four cases are called 
the RG improved gauge action, whose parameters are chosen to be the values suggested 
by approximate renormalization group analyses. Some of these actions are now getting 
widely used, since they realize continuum-like gauge field fluctuations better than the naive 
plaquette action at the same lattice spacing. 

The free gluon propagator is derived in Ref. [loj | : 
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where we employ the Feynman gauge. The matrix A^ v satisfies 



i) A an = for all /i, 



(20) 



(ii) A pu — A up , 



(iii) A pv (k) = A^(-k). 

(iv) V(0) = 1 for n ^ u, 
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and its expression is given by 
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with [i ^ v ^ p ^ r the Lorentz indices. and A 4 are written as 
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In the case of the standard plaquette action, the matrix A^ u is simplified as 

/iplaquette 1 r 

The present calculation requires only the three-gluon vertex which is given in Ref. 
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III. DETERMINATION OF c sw UP TO ONE-LOOP LEVEL 

The first calculation of the clover coefficient up to the one-loop level csw — csw 1 ' '' + 
9 2 csw^ was done by Wohlertjf]], who determined it for the plaquette gauge action to elimi- 
nate the 0(a) contribution in the on-shell quark-quark scattering amplitude. Since the gauge 
propagator is already 0(a) improved, the 0(a) contributions arise only from quark-gluon 
vertex. At tree-level the quark-gluon vertex in Fig. |3] is written as 

A<?>(P, q) = - 9 T A [il, + r (M + M) | _ g r ^L T A £ a, v (p v - q v )a + 0(a 2 ). (34) 

where p and q are incoming and outgoing quark momenta assumed to be much less than 
the cutoff a -1 . We set the Wilson parameter to r = 1. Sandwiching A M (p, q) by the Dirac 
spinor we obtain 

u(q)A^(p,q)u(p) = -gT A u(q){^ + a-c S w {0) )^(p, + q,)}u(p) + O(a 2 ), (35) 

where we use the Gordon identity. We find that csv/ ^ should be one to eliminate the 0(a) 
term. 

To determine the one-loop coefficient csw^ we need six types of diagrams shown in Fig.0] 
The contribution of each diagram to the vertex function is denoted by 

^■\p,q)= E Aj?-">(p, g ) = E f 7^4° (P> <?>*)■ (36) 

i=a,...J i=a,...J J - 7T y Z7T ) 

Here we are concerned with the infrared divergences originating from some types of diagrams. 
Although they are supposed to be canceled out after summing up the contributions of all the 
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diagrams, we need to introduce some infrared regularization in the process of the ca 
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culation. 
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oy the twisted antiperiodic boundary conditions 
7] for this purpose, we instead employ the fictitious gluon 
mass A with the ordinary perturbation theory the infrared divergences are extracted by 
an analytically integrable expression ljp(jp, q, k, A) which has the same infrared behavior as 
If(j?,q,k), 

+ l^^{l^(p,q,k)-9(A 2 -k 2 )iP(j ) , q ,k,X)}\^ o (37) 

with a cut-off A (< n). The Heaviside function 9 is introduced to restrict the domain of 
integration to a hypersphere of radius A, which makes the integral analytically calculable. 
Since we are interested in the 0(g 2 a) contributions, the counter terms lfft(p, q, k, A) can be 
composed of the propagators and vertices, obtained from an expansion of the Feynman rules 
in Sec. illlup to 0(a): 
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D,u(k,X) = ^fp, (43) 
VgixfAh, k 2 , h) = -ipABC {h P (h - k 2 ) T + 2 cycl. perms.} , (44) 

where we consider the massless case. The momentum assignments are depicted in Figs. Q 
and El 

From the Lorentz symmetry and the parity conservation, the off-shell vertex function up 
to 0(p, q) is written as 

+a(P l x + q^Gi + a(p^ - q^Hi + 0(p 2 ,q 2 ,pq)} , (45) 
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where Fi (i = 1, 2, 3), G\ and Hi are dimensionless functions. Sandwiching Ajp(p, q) by the 
on-shell quark states the matrix elements are reduced to be 

u(q)^\p,q)u(p) = -g 3 T A {u(q)-f fl u(p)F 1 + a(p II + q ll )u(q)u(p)G 1 

+a(p» - Q»)u(q)u(p)Hi} , (46) 

where ^w(p) = and u(q)4 — 0. From a view point of the on-shell improvement, the second 
and third terms of the right hand side represent the contributions of the dimension five 
operators, 



+ = (d u ip(x))a llu ip(x) +ip(x)a flv (d u i/}(x)), 
CL = (d v ip{x))o ^{x) -ip(x)a flv (d v tp(x)). 



(47) 
(48) 



Here we should note that the transformation property of CL in terms of charge conjugation 
is different from that of ip(x) , y^ip(x), which means that the last term of eq . (|46j) never appears, 
namely Hi = 0. From the expression ()45|) we can extract the coefficient Gi as 
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It would be instructive to show how the infrared divergence in each diagram cancels out 
after the summation. Let us take the case of the plaquette gauge action as an example. 
Including the constant terms we obtain 
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denotes the contribution of the infrared divergence with the fictitious gluon mass A. The 
integrals are numerically estimated by a mode sum for a periodic box of a size N A with 



N = 64 after transforming the momentum variable through k' = k^ — sin/c^. We choose 
A = 7r for the cut-off. It is found that the tadpole diagram of Fig. 0] (d) gives the dominant 
contribution. The total contribution from infrared divergent terms becomes 

}, (57) 

therefore, the infrared divergences are canceled out in a nontrivial way if and only if the 
tree- level coefficient is properly tuned: csw = 1- Whereas the coefficient of the logarith- 
mic infrared divergence in each diagram is independent of the gauge action, the constant 
terms depend on it. In Table |l] we present the results of csw f° r the various improved 
gauge actions. The value of csw^ 1 "* f° r DBW2 is obtained for the first time. Other re- 
sults are consistent with those obtained by the previous work employing different infrared 
r egular iz at ions []| • 

Here we give a brief description on the mean field improvement of csw- The tadpole 
contribution of Fig. H] is given by 



„ tad _ 2 

csw — 9 



- d 4 k 



(58) 



where fi, v are unsummed and fi ^ v. The numerical values for the various gauge actions 
are listed in Table |l] The mean field improvement is applied as 

c sw = (l + (^C F + JL) s 2 T MF ) (l + S 2 c sw (1) - (JCV + /T MF ) + 0(g A ) 
- i (l + g 2 c sw W ~(\C F + g 2 T MF ) + 0{js% (59) 



u 3 V V3 3N C 

where u = P 1 / 4 is evaluated by Monte Carlo simulation. The derivation of Tmf is given in 
detail in Sec. Ill of Ref . Q- 

The mean-field improved MS coupling g^(fi) at the scale \i is obtained from the lattice 
bare coupling g$ with the use of the following relation: 

IP 22 / 4 \ 

= _ + d g + c p + — log(//a) + Nf [d f - ^ logOxo) J . (60) 

For the improved gauge action one may use an alternative formula 

1 c P + 8c x itt + 16c 2 P2 + 8c 3 P3 

9ms(^) do 
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22 

+d g + (c • c p + 8ci • c R1 + 16c 2 • c R2 + 8c 3 • c m ) + -—r log(/ia) 

l07T 2 

+ JV >( d /-4^ 1 °gW). ( 61 ) 



where 



P = ^TlUplaguette = 1 - Cp#Q + O(^o)> ( 62 ) 
-Rl = 7^TrU rec tangle = 1 ~ C m 0o + O (#o)i (63) 

i?2 = ^Trf/ cW = 1 - c^ 2 + O(^ 4 ), (64) 

-R3 = -Tf U p ara u e l g ram = 1 — Cbz9q + O(^q)) (65) 

and the measured values are employed for P, Rl, R2 and R3. The values of c p , cri, cr 2 and 
c/?3 for various gauge actions are listed in Table XVI of Ref . . 

For later convenience it would be a good idea to parameterize the value of csw^ as a 
function of c\ while keeping c 2 = c 3 = 0. In Fig. EJwe plot the results of csw^ evaluated 
by a mode sum with N = 64, where c\ is chosen from —1.5 to at intervals of 0.02. We 
observe that Csw* 1 ' 1 seems to be divergent as c\ increases. This behavior is well described by 
the rational expression, 

m 0.26849 -0.14193ci- 0.13641c?- 0.07996c? -0.01911c? 
CSW = 1 - 5.08365c, ■ (66) 

where the fitting result is also depicted in Fig. |3J The difference between the actual value 
and the fit is less than 0.1% for —1.5 < c\ < 0. 



IV. CONCLUSION 



In this paper we determine the clover coefficient csw m the massless SW quark action 
up to one-loop order for the various improved gauge actions employing the conventional 
perturbative method with the fictitious gluon mass as an infrared regulator. The validity 
of the method is checked by comparing the results to those previously obtained by the 
Schrodinger functional method: both show a good agreement within error bars. For later 
convenience our results are parametrized in terms of the improvement coefficient c\ of the 
gauge action. An important application of this conventional perturbative method is to 
determine Csw f° r the massive quarks in the heavy quark formulation proposed by the 
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authors, where the relativistic on-shell improvement is extended to the massive case including 
any power corrections of rriQa. Whereas ce and cb receive different mga corrections in 
this formulation, a modification of the present calculational techniques can be done in a 
straightforward manner |l 61] . 
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TABLE I: One-loop coefficient of csw f° r various improved gauge actions. Tadpole contribution of 



Fig. E] (d) is also listed. 



gauge action 


Cl 


C3 


c sw (1) 


_ tad 

csw 


plaquette 








0.26858825(1) 


0.29739454(1) 


Symanzik 


-1/12 





0.19624449(1) 


0.23543879(1) 


Iwasaki 


-0.331 





0.11300591(1) 


0.15988461(1) 


Iwasaki' 


-0.27 


-0.04 


0.12036501(1) 


0.16566349(1) 


Wilson 


-0.252 


-0.17 


0.10983411(1) 


0.15292225(1) 


DBW2 


-1.40686 





0.04243181(1) 


0.08997537(1) 
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FIG. 1: Momentum assignments for the quark- gluon vertices. 
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FIG. 2: Momentum assignment for the three-gluon vertex. 




FIG. 3: Quark-gluon vertex at tree level. 
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(e) (f) 
FIG. 4: Quark-gluon vertex at one-loop level. 
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